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Abstract 

The emergence and spreading of chirality on the early 
Earth is considered by studying a set of reaction- 
diffusion equations based on a polymerization model. 
It is found that effective mixing of the early oceans is 
necessary to reach the present homochiral state. The 
possibility of introducing mass extinctions and mod- 
ifying the emergence rate of life is discussed. Key 
Words: Homochirality - Origin of life - Exobiology. 

Introduction 

The question of the origin of life is closely connected 
with the question of how handedness (chirality) of 
living organisms came about (Avetisov et al. 1991). 
The sugars produced in plants via photosynthesis, for 
example via the rather innocently looking reaction 

6C0 2 + 6H 2 O^C 6 H 12 6 + 6 2 , (1) 

are all right-handed and rotate the polarization 
plane of polarized light in a right-handed sense, i.e. 
the sugar is said to be dextrorotatory. This pre- 
ferred handedness is remarkable because the differ- 
ence in binding energy between right-handed and left- 
handed (levorotatory) sugars is extremely tiny (rela- 
tive energy difference: 10~ 17 ). Once dead, the pre- 
ferred handedness gradually disappears until there 
are equally many right and left handed molecules. 
This process, which speeds up with increasing tem- 
perature, is called racemization and can even be used 
as an approximate dating method (Hare & Mitterer 
1967, Bada et al. 1970). For a simpler and closely 
related reaction that displays a similar behavior we 



refer to the formose reaction discussed by Toxvaerd 
(2005). 

Another type of dextrorotatory sugar with five car- 
bon atoms occurs in the backbone of RNA, for ex- 
ample, making RNA therefore chiral. RNA can have 
autocatalytic properties, which led to the idea of an 
early "RNA world" (Gilbert 1986, Joyce 1991) where 
naked RNA molecules could have catalyzed the poly- 
merization of other RNA molecules. But since the 
RNA molecule consists of right-handed sugars, chi- 
rality selection must have occurred at an even earlier 
stage. 

Two possibilities are commonly discussed. Either 
chirality was imposed by physical factors such as 
magnetic fields, polarized light from a nearby pul- 
sar, or even the tiny energetic preference due to the 
parity-breaking electroweak interaction. In such a 
scenario, homochirality may well have been a prere- 
quisite for the origin of life, allowing for example for 
the assembly of structurally more stable polymers. 
The difficulty here is that such mechanisms would 
hardly explain the complete homochirality observed 
in living matter. The other possibility is that ho- 
mochirality was a consequence of life, and that the 
presently prevailing handedness was selected sponta- 
neously during the assembly of the first polymers. As 
shown in the early paper by Frank (1953), such a pro- 
cess requires the production of left and right handed 
molecules to be mutually antagonistic. He devised 
pairs of ordinary differential equations demonstrating 
how this could possibly be described mathematically. 

The latter scenario is the one we consider also 
in the present paper. Significant progress has been 
made since Frank's early paper. In particular, the 



origin of the mutually antagonistic behavior has been 
identified to be the so-called enantiomeric cross- 
inhibition, a process that spoils further polymeriza- 
tion once a monomer of opposite handedness has been 
attached to an already existing polymer. The evi- 
dence for the occurrence of this process is entirely 
experimental and goes back to the early work of 
Joyce et al. (1984), who found that in non-enzymatic 
template-directed polymerization of RNA strands, 
only a homochiral supply of mononucleotides that are 
complementary to the template can polymerize to a 
typical length of 20 nucleotides. Even a small amount 
of mononucleotides of opposite chirality prevents the 
formation of longer polymers, as is seen in high per- 
formance color chromatograms. Similar experiments 
have subsequently also been carried out by Schmidt 
et al. (1997) and Kozlov et al. (1998), for example. 
In the following we describe in more detail how this 
process can be modeled. 

The polymerization model 

In the polymerization model of Sandars (2003) 
large polymers are generated by joining successive 
monomers into long chains. Monomers can be both 
left, Li, and right, handed, and longer chains 
are formed according to the following set of reaction 
equations: 



L„ + L\ 


2fc S) 


L n +i, 


(2) 




2fcj 


LnRi, 


(3) 


Li + L n Ri 


k S 


L n +iRi, 


(4) 


Ri + L n Ri 


hi 


RiL n Ri, 


(5) 



where k$ and fc/ are the reaction coefficients for 
monomers to a polymer of the same or of the oppo- 
site handedness, respectively. For all four equations 
we also have the complementary reactions obtained 
by exchanging L ^ R. Note that chains are 'spoiled' 
if a monomer of opposite chirality is attached to the 
end of a longer chain, in which case the chain can 
only grow at the other end. This is in essence what 
is meant by enantiomeric cross-inhibition. As a re- 
sults, polymers such as R1L2R1 and can no 



longer grow. As discussed in the introduction, this 
rule is motivated by the experiments of Joyce et al. 
(1984) and others. 

The monomers are initially generated from a sub- 
strate S according to 
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where kc is proportional to the regeneration rate of 
monomers, and Cr and Cl determine the enzymatic 
enhancement of right and left handed monomers. 
The dependence of Cl,r on the existing amount of 
polymer chains is essential to chirality selection as 
then an existing excess of either chirality can be am- 
plified. The exact form of C is not crucial and here 
we follow the choice Ca = J2 n A n , A — L.R (Bran- 
denburg et al. 2005, hereafter referred to as BAHN). 
For alternative prescriptions of Ca we refer to the pa- 
pers by Sandars (2003) and Wattis & Coveney (2005). 
The substrate itself is being replenished by a constant 
source term Q. 

The polymerization process is represented pictori- 
ally in Fig. ^ where one can see how monomers be- 
gin to grow into longer chains and can then be con- 
taminated by a monomer of opposing chirality. The 
crossed out chains represent polymers that can no 
longer grow. 

The feedback mechanism built into the polymeriza- 
tion model leads to an unstable system that, depend- 
ing on the fidelity / of the enzymatic reactions, when 
perturbed from the initial racemic state can reach a 
homochiral state. The prescription for Ca discussed 
above assumes / = 1. If / < 1, there is some 'cross- 
talk' between L and R where / determines the rela- 
tive mixing between the two chiralities. 

The chirality of a particular state is conveniently 
parameterized by a parameter called the enan- 
tiomeric excess: 
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where Ea — J2 n A n with A — L,R. For / larger 
than a critical value (see BAHN for details), the 
racemic state (rj = 0) is unstable with respect to small 
perturbations. Only for / = 1, the end state is fully 
homochiral (rj = ±1). 
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Figure 1: Polymerization process. In addition to a 
number of monomers, L and R, there are several iso- 
tactic dimers, LL and RR, as well as longer polymers. 
Semi-spoiled polymers such as LLLLR and LRRR 
can still polymerize on the unspoiled end. Polymers 
such as LR and RLLR are dead and cannot polymer- 
ize further. 

Homochiralization in space 

The homochiralization process considered here and, 
indeed, in much of the theoretical literature on the 
subject (e.g., Frank 1953, Wei-Min 1982, Goldanskii 
and Kuzmin 1989, Avetisov & Goldanskii 1993, San- 
dars 2003, Saito and Hyuga 2004a) has been uniform 
in space. Relaxing this restriction can, in principle, 
lead to completely new chirality selection mechanisms 
that have no analogue in models without spatial ex- 
tend. An example is the mechanism identified by 
Toxvaerd (2004) using molecular dynamics simula- 
tions. In the present work we simply extend the 
homochiralization process discussed in the previous 
section to allow for the interaction with neighboring 
regions through diffusion and/or advection (such as 
ocean currents). We begin with a state that is ho- 
mochiral everywhere, but that there are minute im- 
balances in space, i.e. infinitesimally small perturba- 
tions. These perturbations grow locally, leading to 
patches with enantiomeric excess of either handed- 
ness. It is then natural to ask how a localized ho- 
mochiral state spreads to a racemic surroundings and, 
moreover, can regions of different chirality coexist? 




Figure 2: Different ways of reaching a homochiral 
state from an initial racemic state 

This question is not restricted to the model consid- 
ered here and needs to be addressed in all models 
where chirality is selected locally by a random pro- 
cess. 

The critical ingredients necessary to answer these 
questions are the presence of different relevant time 
scales. If the appearance of a homochiral region is 
a rare process, say with a period of rrjf e , and the 
global homochiralization time scale, r g i b a i, is short 
compared to that, i.e. r g i b a i <C Tjif c , then one would 
not expect coexisting regions of different chirality. 
In other words, if the emergence of life is a rare 
event, then life forms with different handedness did 
not probably coexist. This is the case modeled by 
Frank (1953) and many others after him. In the other 
extreme case, homochiral regions appear frequently 
compared to the speed at which they can dominate 
the early Earth, leading to coexisting regions of op- 
posite chirality (see Fig. Even though here we 
are mainly concerned with chirality, such arguments 
can be applied more generally to the emergence and 
spreading of life as well. For a discussion of the pos- 
sibility of finding a second sample of life on Earth see 
the recent paper by Davies & Lineweaver (2005). 

The spreading of chirality and coexistence of re- 
gions of different chirality was studied in Branden- 
burg and Multamaki (2004), hereafter referred to 
as BM. There a reduced polymerization model with 
spatial extent was considered analytically and by 
numerical simulations. The importance of spatial 
extent has already been emphasized by Saito and 
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Hyuga (2004b) who generalized the model of Saito 
and Hyuga (2004a) by using a Monte-Carlo method 
as it is used in percolation studies. 

The behavior of the reduced model can be under- 
stood qualitatively: an initial racemic mixture with 
small local perturbations quickly relaxes locally to 
regions of different chirality. Hence, in this work the 
emergence of life was considered to be a frequent and 
rapid event. The oppositely handed regions then be- 
gin to spread spatially into any possibly remaining 
racemic regions by front propagation until they come 
in contact with a region of opposite chirality. In 1 + 1 
dimensions (one space and one time dimension) there 
is no further evolution and one can view the process 
as spontaneous symmetry breaking leading to sta- 
ble, non-propagating domain walls. In more than one 
spatial dimension, the homochiralization process pro- 
gresses further. If no advection is present, only diffu- 
sion can drive the homochiralization process. Analyt- 
ical arguments can be utilized to show that a bubble 
surrounded by a region of opposite chirality tends 
to shrink. Since the equations are local, the local 
curvature of the front is the deciding property in de- 
termining which way a front will move. An example 
of this process is shown in Fig. |3 where we show a 
2 + 1 dimensional box with regions of different chiral- 
ity. The initial state in the simulation was a racemic 
state with small fluctuations so that one quickly ar- 
rives at a state where there are many left and right 
handed regions in the box. Here we show the evolu- 
tion of these homochiral regions. From the figure it 
is easy to see that the curvature of the interface in- 
deed determines the evolution of a region: e.g. small 
bubbles shrink and later disappear. 

Again, time scales play a crucial role in the global 
homochiralization process. Using microphysically 
motivated values, we find that the local relaxation to 
a homochiral state is rapid compared to that of the 
propagation of fronts into racemic surroundings. In 
turn, this is a much faster process than the final stage 
in global homochiralization when the regions of dif- 
ferent chirality slowly disappear as described above. 
Using values of molecular diffusion, one actually finds 
that this time scale on Earth is much too long to ex- 
plain the observed homochiral final state. 

In order to study the process further, we have also 
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Figure 3: The spreading of homochiral regions in 
2 + 1 dimensions. The dark (light) regions correspond 
to left (right) handed regions. Time is given in global 
diffusion times, r tU rb = L 2 / n. (Adapted from BM.) 



considered the effects of advection on homochiraliza- 
tion. An example in 2 + 1 dimensions is shown in 
Fig. The mixing process is now greatly enhanced 
and one finds that the global homochiral state is 
reached much more rapidly. The actual time scale 
depends strongly on the strength of the flow, e.g. for 
root mean square flow of 1 cm/s, the time scale of 
global homochiralization is of the order of 30 years. 
It is then clear that if the model considered in BM 
captures the relevant features of chirality selection on 
the early Earth, the effective mixing of early oceans 
is vital. As the mixing is affected by many factors 
such as the existence of continents and salinity of sea 
water, it is not difficult to speculate that in some se- 
cluded parts of the early oceans, life forms of different 
chirality could have coexisted. 

So far we have considered the spreading in two hori- 
zontal directions, which is relevant to the Earth's sur- 
face where the vertical dynamics may be eliminated 
by vertical averaging. However, BM also considered 
three-dimensional models and found that the speed of 
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Figure 5: The spreading of homochiral regions in 3 + 1 dimensions. The dark (light) regions correspond to 
left (right) handed regions. The time t = 50 corresponds here to t/rturb ~ 2. This is also the last time shown 
in Fig. 0J By the time t — 100 corresponds here to t/rturb ~ 4, the right handed life form went extinct. 
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Figure 4: The spreading of homochiral regions in 
2 + 1 dimensions. The dark (light) regions corre- 
spond to left (right) handed regions. Time is given 
in dynamical times, r tur b = i/u rms . (Adapted from 
BM.) 



homochiralization is somewhat enhanced (by about 
30%) relative to the two-dimensional case. Further- 
more, in both two and three dimensions they found 
that the generalized enantiomeric excess, 
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where angular brackets denote averages in space, 
grows approximately linearly in time: 



^ = (12. ..18) x J\f 
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where Af is the number of topologically disconnected 
domains, X/yQlcs is the non-dimensional growth 
rate in the spatially uniform case, L is the spatial 
extent of the domain, and K is the diffusion coeffi- 
cient. 

In the turbulent case, k can simply be replaced 
by a turbulent coefficient Kturb, for which BM found 
from their simulation the rough estimate Kturb ~ 
(0.08...0.15)u rms £, where u rms is the root mean square 
velocity of the turbulence and £ is the typical length 
scale of the turbulent eddies. It should be empha- 
sized, however, that such a description of turbulent 
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mixing is rather rudimentary and ignores, for exam- 
ple, the fact that the shape of left and right handed 
domains is quite different with and without turbu- 
lence: it is more round in the case of molecular diffu- 
sion and more elongated in the turbulent case. The 
latter is a feature of the turbulence to mix fields by 
stretching and folding. A snapshot of a fully three- 
dimensional simulation 1 showing the concentration of 
left handed building blocks is presented in Fig. |SJ 

Emergence and spreading of life 
on the early Earth 

We have made a number of underlying assumptions 
in studying the spreading of homochiral domains on 
the early Earth. Firstly, we have assumed that the 
initial state is globally racemic. This implies that the 
conditions were life could have emerged are uniform 
and favorable everywhere and, furthermore, that a 
homochiral region can spread and fill the whole do- 
main. Secondly, we have assumed that the emergence 
of life is a common event and that it is possible ev- 
erywhere on earth, and that there are no mass ex- 
tinctions. These assumptions may not necessarily be 
realistic, and one should consider what new effects 
can possibly arise from relaxing these assumptions. 

Considering the possibility of mass extinctions and 
that the emergence of life can be a rare event, it is 
clear that a new time scale, r ex t, is inserted into the 
process. Assuming that mass extinctions are local 
events (i.e. not all life is destroyed in which case the 
process can start again), the relative lengths of the 
different time scales determine the qualitative dy- 
namics of the system. If r g i b a i *C ruf e <C T cxt , life 
appears somewhere and takes over the whole system 
before life of opposite chirality can emerge. Any later 
mass extinction events will be insignificant as the 
dominant chirality will quickly win over the achiral 
region. If ri ifc <C T glob ai -C r GXt , we have the situa- 
tion considered earlier in this paper, where regions of 
different chirality compete. If mass extinctions are 

1 Animations of both turbulent and 

purely diffusive solutions can be found at 
http : //www . nordita . dk/~brandenb/movies/chiral 



common, r oxt < (rnf c , T g i oba i), and so any emergent 
life will quickly be wiped out - at least locally. Again, 
if ^global "C Tiif e , it is unlikely that life will re-emerge 
spontaneously in the affected areas, and these areas 
will more likely be re-populated by the spreading of 
the homochiral regions surrounding the now racemic 
area. Thus, in this case homochirality is preserved. 
On the other hand, if T\if e <C Tgi ODa i, the possibility of 
mass extinctions allows for new life forms to emerge 
which, in turn, may prolong the time during which 
life forms of opposite chirality can have coexisted. 

It is not immediately obvious how mass extinc- 
tions should be modeled within the framework of the 
BAHN model. It seems plausible that a mass extinc- 
tion due to an impact, for example, can be modeled 
by removing all existing polymers. The re-emergence 
of life depends on the source term Q of the substrate 
and on the kc parameter that controls the rate at 
which new monomers are created from the substrate. 
By locally suppressing Q or k c as a consequence of 
the mass extinction, one can model the reduced hab- 
itability by slowing down the generation either of the 
substrate or of new monomers from the source. Al- 
ternatively, one could decrease the polymerization ef- 
ficiency by decreasing ks and ki on the grounds that 
these parameters are likely to depend on external fac- 
tors such as temperature, salinity, and acidity. 

Allowing for the additional dependencies discussed 
above introduces certainly a lot of additional uncer- 
tainties, but it also removes some of the sensitivity to 
initial conditions. If the conditions for chirality (or 
life) are not satisfied everywhere, the spreading of 
chirality will stop and mass extinctions can eradicate 
the emergent life forms. Hence, there are various ap- 
proaches to modeling the emergence of life and mass 
extinctions in the BAHN model. Either one can de- 
plete the source of the substrate, or one can diminish 
the conversion rates for monomer formation and/or 
for polymerization, in which case the substrate might 
still be present everywhere. 

Clearly, these considerations are relevant in a more 
realistic model and, furthermore, can bring interest- 
ing new dynamics into the system. As such they are 
worth considering in more detail and future work will 
hopefully bring us answers to these questions. De- 
tailed numerical models of this process may elucidate 
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further the possible outcomes. We plan to adapt the 
Pencil Code 2 , which was also used in our earlier 
paper, to include the effects of a dynamical evolu- 
tion of the rate at which new left and right handed 
monomers can be regenerated. 

Conclusions 

We have studied the spatial evolution of chirality in 
an initially racemic system. By utilizing analytical 
and numerical methods, we have shown how a glob- 
ally homochiral state can be reached. The associated 
time scale depends strongly on the details of the sys- 
tem, but molecular diffusion is too slow to homochi- 
ralize the early Earth and turbulent flows are neces- 
sary. Effective mixing of the early oceans drastically 
reduces the required time and is hence a vital ingredi- 
ent in understanding the homochiralization process. 

We also consider qualitatively how a more infre- 
quent emergence of homochiral regions and the ad- 
dition of mass extinctions would affect the evolution 
of the system. Adding such (possibly more realistic) 
processes enriches the dynamics and offers new in- 
teresting avenues of research to be explored in the 
future. 
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